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^ ! ABSTRACT 

CN) ' Context. We present observations of rotational lines of H2S, SO and CS performed in comet C/1995 Ol (Hale-Bopp) in March 1997 

with the Plateau de Bure interferometer of Institut de Radioastronomie Millimetrique (IRAM). 

Aims. The observations provide informations on the spatial and velocity distributions of these molecules. They can be used to constrain 
their photodissociation rate and their origin in the coma, i.e. nucleus or parent source. 

Methods. We use a radiative transfer code which allows us to compute synthetic line profiles and interferometric maps, to be compared 
to the observations. 

\Zi • Results. Both single-dish spectra and interferometric spectral maps show evidence for a day/night asymmetry in the outgassing. From 

' the analysis of the spectral maps, including the astrometry, we show that SO and CS present in addition a jet-like structure that may 

' be the gaseous counterpart of the dust high-latitude jet observed in optical images. A CS rotating jet is also observed. Using the 

astrometry provided by continuum radio maps obtained in parallel at IRAM, we conclude that there is no need to invoke the existence 
of nongravitational forces acting on this comet, and provide an updated orbit. The radial extension of H2S is found to be consistent 
with direct release of H2S from the nucleus. SO displays an extended radial distribution. Assuming that SO2 is the parent of SO, the 
^ ' photodissociation rate of SO is measured to be 1 .5 x lO"** s"' at 1 AU from the Sun. This is lower than most laboratory-based estimates 

04 , and may suggest that SO is not solely produced by SO2 photolysis. From the observations of 7(2-1) and J(5^) CS lines, we deduce 

0^ ■ a CS photodissociation rate of 1 - 5 x 10"^ s"'. The photodissociation rate of CS2, the likely parent of CS, cannot be constrained due 

' to insufficient angular resolution, but our data are consistent with published values. 

CO ' Conclusions. These observations illustrate the cometary science that will be performed with the future ALMA interferometer. 

0^ ' Key words. Comet: individual: C/1995 Ol (Hale-Bopp) - Radio lines: solar system - Techniques: interferometric 

o ■ 

. . 1. Introduction iFeldman et al.ll2004l) and provided measurements of their pho- 

todissociative lifetime and the lifetime of their parents. Spatial 
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Spectroscopic investigations of comets in the microwave wave- ^.j^ high angular resolution has also been under- 

lengths domain began in the 80 s when sensitive instrumenta- ^^^^^ ^ molecules (CO, HCN, H2O) observed at infrared 



^ tion was made available. This wavelength range proved to be a wavelengihs"in"comet"s C/1996 B2 (Hykutake )"and/or "c/"r995 

- - powerful tool to identify and measure the relative abundances of pi (Hale-Bopp) (e.g., DiSanti et aL .2001: Magee-Sauer et all 

parent molecules released from comet nuclei from the observa- 1^9991) 
tions of their rotatio nal lines (see, e.g., the review of Bockelee- 
Morvan et al. 12004 . However, little information was obtained 

on their spatial distribution in the inner coma because of the Mapping of rotational emission lines at high angular reso- 
poor angular resolution (typically 10" at most) offered by sin- lotion was only attempted in a few comets, mosdy in comets 
gle dish (diffraction-limited) measurements. The study of the Hyakutake and Hale-Bopp. Because of their exceptional bright- 
spatial disti-ibution of parent molcules is of sti'ong interest as "^ss, interferometric imaging of molecular lines could be suc- 
inferences concerning the outgassing processes at the nucleus cessfuUy performed (e.g., | Wright et al. || 1998HBlake et alj|1999| ; 
surface can be made. In addition, the study of the brightness ra- IWoodney et al.|2002|) . We present here observations of three sul- 
dial distribution of molecular species can provide information fur bearing species undertaken in comet Hale-Bopp in March 
on their origin, i.e. indicate whether they are indeed released 1997 using the IRAM Plateau de Bure interferometer at an- 
by the nucleus or produced by some source in the coma. Such S^^^^ resolution of 1.6 to 3.5". H2S and SO have so far 
investigation has been performed for decades on radicals and only been detected m comets in the rnillimetre or submil- 
atoms observed at visible or UV wavelengths (see the review of limeti-e domains (|Bockelee-Morvan et al.l |2004f). Our observa- 
tions provide the first study of their spatial distribution. CS 



Send ofprint requests to: LBoissiei has been detect ed at h igh angular resolution in the UV in 

* Based on observations carried out with the IRAM Plateau de Bure several comets dPeldman et al. 2004). Its spatial distribution 

Interferometer. IRAM is supported by INSU/CNRS (France), MPG was studied in the microwave domain in comet Hale-Bopp us- 

(Germany) and IGN (Spain). ing the Berkeley Illinois Maryland Association (BIMA) array 
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dSnvder et al.l200lT) . Our observations of two CS rotational lines 
usefully complement these studies. 

Section|2]presents the data set. Interferometric images, spec- 
tral maps and autocorrelation (i.e., single-dish) spectra of H2S, 
SO and CS lines are analysed and compared to investigate jet- 
like structures and temporal variability linked to nucleus rota- 
tion. In Sect. [3] we show that a well-determined gravitational 
solution for comet Hale-Bopp orbit can be deduced from the as- 
trometry provided by the Plateau de Bure continuum observa- 
tions. Section|5]focuses on the study of the radial distribution of 
the molecules, the model being described in Sect.E] A summary 
follows in Sect. |6] 



2. Observations 

2.1. Description 

Comet C/1995 Ol (Hale-Bopp) was observed in 1997 at the 
IRA M Plateau de Bure Inte rferometer situated in the french 
Alps dGuilloteau et al.L Il992h . The observing campaign lasted 
2 weeks between March 6 and March 22, at the time of the 
closest approach of the comet to the Earth (perigee was on 
March 22). The programme included interferometric mapping 
and single dish measurements in spectroscopic mode for sev- 
eral molecular spe cies (CO, HCN, HNC, CS, H2S, SO2, SO 
CH3OH, H2 CO; Wink et al. 1999;'Bockelee-Morvan et al. 2000t 
lHenrvetalJl2002; Henry ,2003; Bockelee-Morvan et al. 2005^ , 
as w ell as observations of t he dust and nucleus continuum emis- 
sion (lAltenhoffetalJll999h . 

At the IRAM Plateau de Bure interferometer, observations 
with 1.3 and 3 mm receivers were conducted at the same time. 
The CS 7(2-1) and 7(5-4) lines (at 98.0 and 244.9 GHz, respec- 
tively) were observed in parallel on March 12. On March 13, the 
22o-2ii line of H2S (at 216.7 GHz) and A^y(56^5) line of SO 
(at 219.9 GHz) were observed simultaneously in the same back- 
end, taking benefit of double-sideband receivers. A log of the 
observations is presented in Table[T] At this time, the instrument 
counted five 15-m antennas set in CI compact configuration so 
that baseline lengths projected onto the plane of the sky were be- 
tween 22 m and 147 m. The comet was at a geocentric distance 
A = 1.35 - 1.36 AU and heliocentric distance rn = 0.97 - 0.98 
AU. 

The comet was tracked using orbital elements computed by 
D.K. Yeomans (JPL ephemeris solution 55) and the ephemeris 
program of the Institut de mecanique celeste et de calcul des 
ephemerides (IMCCE, Observatoire de Paris) that takes into 
account planetary perturbations. First interferometric maps ob- 
tained on March 9 showed that the intensities (of both continuum 
and line emissions) peaked 5-6" northward with respect to the 
expected position. The ephemeris was corrected by +6" in dec- 
lination (Dec) for the observations performed on March 12 and 
13. It turned out that the position of the brightness peaks (con- 
tinuum and lines) was offset by less than ~2" with respect to 
the pointed position. Single dish measurements are not affected 
by this pointing offset given the size of the primary beam at the 
observed frequencies (Half Power Beam Width, HPB W, of 22 to 
50", Table[B. 

The observing cycle was typically: pointing, focussing, 4 
min of cross-correlation on the calibrators, 2 min of autocor- 
relation (ON-OFF measurements, 1 min on source) and 51 
scans (1 min each) of cross-correlation on the comet. At the 
end of some cycles, additional autocorrelation spectra were ob- 
tained. Hereafter, autocorrelation observations will be refered to 



as "ON-OFF" observations. The OFF position in ON-OFF (po- 
sition switching) observations was at 5' from the ON position. 

In 1997, the correlator consisted of 6 independent units, ad- 
justable in resolution and bandwidth. On March 12, two units 
were used for observing CS, the other ones were used for 
CH3OH lines near 241 GHz. On March 1 3, four units were use d 
for measuring the continuum emission (lAltenhoffetal.1119991) . 
and one unit c overed the 7(1-0) HNC line o bserved with the 3 
mm receivers (iBockelee-Morvan et al.ll2005l) . Another one cov- 
ered both the 22o-2ii H2S and NjiS^-^s) SO lines present in 
lower and upper receiver sidebands. The units used for CS, H2S 
and SO lines were divided into 256 channels with 78 kHz separa- 
tion, providing an effective velocity resolution of ~ 0.3 1 km s"' 
forCS 7(2-1), ~ 0.12 km s~' forCS 7(5-4) and ~ 0.13 km s"' 
for H2S and SO lines (the resolution is 1.3 times the channel 
spacing). 

The frequency calibration was checked by observing the ref- 
erence source W30H. Interferometric observations were cali- 
brated using the standard IRAM software CLIC. The amplitude 
and phase calibrator was BL Lac (2200+420) which flux density 
was determined observing MWC349. In ON-OFF mode, the in- 
tensity scale is the antenna temperature Ta*. The conversion fac- 
tor from Ta* to flux density and the HPBW of the primary beam 
at the frequencies of CS, H2S and SO lines are given in Table[T] 
The uncertainty in flux calibration is about 15%, both in inter- 
ferometric and ON-OFF mode. 

2.2. ON-OFF spectra 

The ON-OFF spectra obtained independently by the five anten- 
nas were co-added. All lines were detected with a signal-to-noise 
ratio > 10 in single scans, so that the evolution of the line char- 
acteristics along the observation period could be studied to in- 
vestigate possible changes linked to nucleus rotation. Since the 
rotatio n period of the nuc leus of comet Hale-Bopp is about 1 1 .35 
h (e.g.. lJorda et al.lll999l) . the nine H2S and SO spectra obtained 
between 4hl9 and 14h37 UT (Table[TJ spread over almost a full 
rotation period. Similarly, the observations of CS 7(5^) cover 
80% of the nucleus period. Technical problems prevented us 
from obtaining as much time coverage for CS 7(2-1). Spectra 
as a function of time are shown in Fig. [TJ[2l 

As shown in Figs. 1-2, no significant temporal changes that 
could be related to nucleus rotation are observed for the line ar- 
eas. We also plot in Figs. 1-2 the temporal evolution of the line 
centroids in nucleus velocity frame, hereafter referred to as line 
velocity shifts or velocity offsets (see footnote (c) of Table 1). 
For all molecules, velocity offsets are always negative, taking 
into account error bars. Given the phase angle on March 12-13 
(~ 45°), this is consistent with a Sun/anti-Sun asymmetry in the 
distribution of the molecules, with more molecules towards the 
Sun. The velocity offsets of the H2S (22o-2i 1), SO A^y(56-45) and 
CS 7(2-1) lines present some fluctuations but no clear evidence 
for a 11.35 h periodic modulation (Figs. [T]-|2]l. In contrast, the 
velocity offset of the CS 7(5^) line, measured with much bet- 
ter accuracy (less than 0.01 km s ' rms), presents a clear trend 
of variation with rotation phase, with a minimum on March 12 
10hl5 UT at -0.125 km s"', and maximum at -0.03 km s^' 
observed 4 h earlier 

These measurements have to be compared to the strong mod- 
ulation (~ 0.3 km s"' total amplitude) observed on March 1 1 for 
the CO 7(2-1) and 7(1-0) velocity offsets (the lowest velocity 
offset of the CO 7(2-1) line was near 8h30 UT). This modula- 
tion has been interpreted as caused by a rotating CO jet emanat- 
ing from a source near the nucleus equator that was active day 
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Line 


Date 


T TT 
U i 




A 


ivioae 


^ 1 1 A 


Beam 


Line area 


AV 
















\ ) 


Jy km s 


km s 


CS 7(2-1) 


March 12 


8h47-15h30 


0.98 


1.36 


cross 


23.5 


4.39 X 2.82 


1.47 ± 0.09 


-0.19 ± 0.02 












auto 


9 1 1 

ill . 1 


AO 


LJ.A It U.Z, 




CS 7(5-4) 


March 12 


6h31-15h30 


0.98 


1.36 


cross 


63.1 


1.92 X 1.28 


5.4 ±0.8 


-0.10 ± 0.02 


244.9355565 




6h24-14h03 






auto 


30.4 


19.7 


213.2 ±0.6 


-0.079 ± 0.003 


H2S 22o~2[] 


March 13 


4h25-15h27 


0.97 


1.35 


cross 


37.6 


1.99 X 1.53 


3.51 ± 0.15 


-0.05 ± 0.01 


216.7104350 




4hl9-14h37 






auto 


27.0 


22.15 


18.9 ±0.3 


-0.06 ± 0.02 


SO A'y(56~45) 


March 13 


4h25-15h27 


0.97 


1.35 


cross 


37.6 


1.93 X 1.5 


0.44 ± 0.04 


-0.05 ± 0.09 


219.9494420 




4hl9-14h37 






auto 


27.0 


21.9 


16.9 ± 0.3 


-0.13 ±0.02 



Conversion factor averaged over the six antennas. Forward and beam efficiencies are (Feff, B^ff) = (0.93, 0.83), (0.89, 0.62), (0.89, 0.55) at 
98, 220 and 245 GHz, respectively. 

Intensity integrated over velocity measured in ON-OFF spectra (autocorrelation mode), and in the central pixel of the interferometric maps 
(cross-correlation mode). The given uncertainties do not include uncertainties in the absolute calibration (~ 15% in both modes). 
Line velocity shift defined as Av = y^'' , where T, and v, are the intensity and velocity of channel respectively 



and night jHenrv et alj|2002l: lHenrvll2003l; Bockelee-Morvan et 
al. in preparation: \Boissier et al.ll2006 '). as in mid-March the nu- 
cleus rotation axis was near the plane of the sky. The velocity 
offset of the CO lines took negative to positive values because 
the velocity vector of the CO-jet with respect to Earth rotated 
during the course of the observations. Comparing CS and CO 
lines evolution, one can conclude that their evolution are not in 
phase, and thus that the CO and CS sources might be different. 

The average spectra obtained on March 12 and 13 are shown 
in Fig. [3] Their characteristics are summarized in Table [T] 
Comparing the shapes of the lines and their Doppler shift, we 
can notice that the SO line is more asymmetric than the H2S 
line. This will be further discussed in Sect. 12.31 



2.3. Interferometric observations 

2.3.1 . Principle of interferometric imaging 

An interferometer measures the 2D Fourier Transform (FT) of 
the brigthness distribution on the sky. At a time t a pair of anten- 
nas provides a point of the FT (a complex number called visibil- 
ity) at the coordinates (u, v) in the Fourier Plane (also called uv- 
plane), corresponding to the projection on the sky of the baseline 
(separation between antennas). The baseline projection evolves 
due to Earth rotation so that a single pair of antennas probes sev- 
eral points of the FT during the observations. At the end of the 
observing run, the visibilities are stored in a Mv-table. The final 
observed sample in the Fourier plane is called uv-coverage. 

Interferometric maps are obtained by computing the inverse 
Fourier Transform of the measured signal. However, because in 
practice it is impossible to probe the whole Fourier plane, inter- 
polation is required. The asymmetry in the uv-coverage results 
in an elliptic clean beam. The clean beam (hereafter referred to 
as interferometric beam) is the central part of the FT of the uv- 
coverage. Direct analysis of the visibilities allows us to avoid the 
uncertainties due to the interpolation. 



2.3.2. CS, H2S and SO interferometric maps 

The interferometric maps, obtained by co-adding the velocity 
channels where line emission is present, are shown in Fig.|4] The 
original images have been shifted so that the maximum bright- 



ness peaks at the centre of the maps. The main characteristics of 
these maps (intensity at maximum brightness, absolute and rel- 
ative coordinates of the peak brightness (O) with respect to cur- 
rent ephemeris (C) and with respect to the position of maximum 
continuum emission (Cont) measured simultaneously) are sum- 
marized in Tables [1] and |2] and Fig.|5] Spectral maps are shown 
in Figs. |6H8] The following remarks can be done : 

- The SO spatial distribution is more extended than that of 
H2S. The line flux in the interferometric beam is indeed eight 
times lower for SO (Table [TJ, while SO and H2S lines have 
comparable intensities in the primary beam, as measured in 
autocoiTelation mode (Fig. |3] and Table [TJ. Radial distribu- 
tions will be studied in detail in Sect. |5] 

- The points of maximum line brightness do not coincide with 
the nucleus position provided by the ephemeris (C), found 
to be 1.5 to 3" southward (Fig.|3. Similar (O-C) offsets (by 
typically 1" in RA and 3" in Dec) are observed for other 
lines (CO , HCN, HNC, H2CO CH3OH) mapped at Plateau 
de Bure (Altenhoff et al.l Il999 l, Bockelee-Morvan et al. in 
preparation). Continuum maps obtained simultaneously at 
Plateau de Bure also show a significant trend to peak north- 
ward o f the ephemeris positi on, with (Cont-C) offsets in Dec 
~+3" (lAItenhoffetal.l(l999l) . 

- The positions of H2S and CO maximum emission match 
the continuum position. No continuum data were obtained 
on March 12 at the date of CH3OH observations, but al- 
most identical (Cont-C) offsets were measured on March 9, 
11 and 13 (Fig.|5]l. Assuming that this (Cont-C) offset ap- 
plies on March 12, the methanol emission peak also matches 
the continuum position. In contrast, the point of maximum 
brightness of other molecules (namely HCN, H2CO, CS, SO 
and HNC) lie at intermediate positions between the contin- 
uum and the ephemeris. 

- In the spectral maps (Figs. ISHSJ, the lines exhibit larger 
blueshifts at northward positions with respect to the centre 
of brightness in the maps. As already noted, blueshifted lines 
are expected for preferential outgassing towards the Sun. 
Taking into account the Earth-Comet-Sun geometry, stronger 
blueshifts towards northern positions is also consistent with 
a day/night asymmetry. Indeed, though the Sun is towards 
South in the maps (PA ~ 160°), a larger number of molecules 
released from the nucleus illuminated hemisphere are ob- 
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Table 2. CS, H2S and SO maximum line brightness positions. 



Line 


Date 


UT" 


RA'' 


Dec'' 


(0-cont)'^ 


(O-C)'' 












((5RA,5Dec) 


((5RA,(5Dec) 






h:min:s 


h:min:s 


Q .1 .tt 


("/') 


(",") 


CS i(2-l) 


March 12 


15:00:00 


22:41:28.643 


41:36:31.11 


e 


(+0.42,+2.13) 


CS 7(5-4) 


March 12 


7:00:00 


22:38:43.424 


41:24:11.52 


e 


(+0.37,+2.48) 


H2S 


March 13 


8:00:00 


22:47:25.371 


42:02:02.88 


(+0.06,4-0.05) 


(+0.92,+2.87) 


SO 


March 13 


5:00:00 


22:46:21.868 


41:57:35.54 


(-0.03,-1.34) 


(+0.83,+ 1.48) 



^ UT time at which coordinates are given. 
Apparent geocentric coordinates of the brightness peak measured on interferometric maps ("O"). Gaussian fits in the Mv-plane provide similar 
absolute positions. 

° Offset in RA and Dec between molecules brightness peak position ("O") and 1.3-mm dust continuum peak position ("cont". IXltenholf et al.l 
Continuum observations were not performed on March 12. 

^ Offset in RA and Dec between molecules brightness peak position ("O") and the comet expected position ("C") based on astrometric mea- 
surements from 27 April 1993 to 4 August 2005 (JPL ephemeris solution 220). 
Continuum observations were not performed on March 12 CAltenhoff et al...l999,) . 



served to outflow northward rather than southward because 
the phase angle is ~ 45°. This was checked by synthesiz- 
ing spectral maps using the model detailed in Sect. |4] and a 
Sun/anti-Sun aymmetry in the distribution of the molecules 
in the coma. 

- There are however significant differences between H2S and 
SO and CS spectral maps. The H2S spectra at northern po- 
sitions are less asymmetric than those of SO and CS at the 
same positions. In addition, the CS 7(5-4) spectra at centre 
and southern positions exhibit excess emission at Doppler 
velocities around km s"'. At these velocities, the emitting 
molecules are moving in (or near) the plane of the sky. 

The molecules can be then grouped in two categories: those 
whose brightness centres coincide with the continuum (H2S, 
CO, CH3OH) and whose spectral maps (averages over about one 
rotation period) show mainly evidence for a global day/night 
asymmetry. Molecules of the second group (CS, SO, but also 
HCN, HNC and H2CO) have their brightness centres offset to- 
wards South with respect to the continuum and may present, in 
addition to a diurnal asymmetry, a jet-like structure outflowing 
towards South in the plane of the sky, as observed for CS on 
its 7(5-4) spectral maps. Comet Hale-Bopp exhibited in March 
1997 a strong du st jet pointing southwest in visible images 
jjorda et alJI 19991) and interpreted as originating from high lat- 
itude regions on the nucleus surface (the rotation axis was at 
PA ~ 210° and aspect angle ~ 70° in rnid-M arch 1997, i.e. al- 
most in the plane of the skv. lJorda et al.l ll999). We suggest that 
molecules of the second group present a gaseous counterpart of 
this dust jet. 

On the basis of synthesized maps obtained by describing 
the coma with a combination of conical structures of various 
opening angles and, for the continuum, a point-source mimick- 
ing nucleus thermal emission, we conclude that: 1) the strong 
offset (Cont-C) ~ 3" in Dec cannot be explained by asymme- 
tries in the dust coma, assuming a 50% nucleus contribution to 
the thermal flux as derived by Altenhoff et al. ( 1999); the cen- 
tre of brightness of the continuum maps should be very close 
(within 0.3") to the actual position of the nucleus; 2) Sun/anti- 
Sun asymmetries in the outgassing that are compatible with ob- 
served blueshifts cannot displace the centre of molecular line 
brightnesses by more than a fraction of arcsec, thereby explain- 
ing why molecules of the first group have their centre of bright- 
ness almost coinciding with the continuum peak (i.e. near nu- 



cleus position); 3) the (O-Cont) offsets in Dec observed for the 
molecules of the second group can be explained by enhanced 
production in a high-latitude jet-like structure; for CS, a 0.8" 
off'set is obtained with a model including a wide (~ 50°) coni- 
cal molecular jet oriented in the plane of the sky and comprising 
about 25% of the total production of CS (Fig. |9]i; for SO, the 
1.4" offset is obtained with a wider jet (~ 80°) containing about 
40% of the total production. Such a static wide jet could cor- 
respond to the temporal average of a narrower rotating jet; the 
simulations show that a model combining a diurnal asymmetry 
and a high-latitude jet reproduces qualitatively the line shapes of 
CS and SO ON-OFF spectra and spectral maps as well as the 
(O-Cont) offsets (Fig.|9]l. 4) the position of the centre of bright- 
ness in time-integrated maps is not significantly affected by low- 
latitude rotating gaseous jets, such as the one observed for CO 
(Sect.|221l. 

Fitting in detail the H2S, CS and SO interferometric maps 
is challenging, given the complexity of Hale-Bopp coma and 
lack of temporal coverage. This study is beyond the scope of the 
present paper In Sect.|5] we will focuss on the radial distribution 
of these molecules. 



3. Astrometry: orbital implications 

It is at this point useful to discuss the Plateau de Bure obser- 
vations in the context of the orbit determination of the comet. 
Whether it is or is not necessary to allow for effects of other 
than gravitational forces o n the Hale-Bopp orbit has long been 
a rather contentious issue (lMarsdenl[T999E . On the one hand, it 
was difficult to believe that, active though the comet was, the rel- 
ative mass loss due to outgassing from a nucleus as large as Hale- 
Bopp was estirnated to be (from 45 to perha ps 70 km diameter, 
lAltenhoff et al.ll999llWeaver & Lamvll999i) would be sufficient 
to produce a detectable effect on the orbit. On the other hand, as 
the timespan of the astrometric observations increased well be- 
yond the comet's perihelion passage, the manner in which the 
orbital solutions satisfied the data was becoming less and less 
satisfactory, even when special weighting schemes were intro- 
duced to stress some observations over others. Eventually, when 
the observations extended some ten months after perihelion pas- 
sage, there was no other possibility but to incorporate nongrav- 
itational parameters into the solution, including in particular a 
radial force component outward from the Sun that had a formal 
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mean error of only some three percent. In broad terms, later or- 
bital solutions, whether computed at the Minor Planet Center or 
at the Jet Propulsi on Lab oratory (e.g., MPC 54799 and JPL 220 
orbit used in Sect 12.3.21 which both include observations made 
into Aug. 2005), seemingly confirmed this finding. 

If there were a significant departure between the centre of 
mass and the centre of brightness when the comet was most ac- 
tive, it makes sense to attempt entirely gravitational orbital solu- 
tions that exclude consideration of the optical astrometry during 
that time period. Although the last observations were made more 
than 12 years after a fortuitous prediscovery observation, some 
19 percent of the total of 3565 observations in the Minor Planet 
Center's files were made during the 50 days immediately preced- 
ing the comet's perihelion passage. The exclusion of these obser- 
vations permits a moderately satisfactory gravitational fit to the 
remainder. The (O-C) residuals in Dec for the omitted times- 
pan are systematically negative to an average extent of about 
3". Only 2 percent of the observations were made during the 
two months after perihelion passage, but this negative trend also 
shows during this period of included observations, after which 
there was a six-week gap in the observational record because 
of the comet's small elongation from the Sun. There is likewise 
a negative trend during the last month of included observations 
prior to perihelion passage. 

Orbit I in Table[3]was therefore computed from 2693 obser- 
vations covering the time periods 1995 July 24-1997 Jan. 10 and 
1997 July 14-2005 Aug. 4. Observations with a residual in RA 
or Dec or both of 2.0" or larger were also excluded from this 
gravitational solution. The mean residual of the included data 
is 0.86", and there are no obvious systematic trends. The sin- 
gle prediscovery observation on 1993 Apr. 27 has an (O-C) Dec 
residual of -2.0", and the negative trend during the omitted half- 
year is much as before. 

The ten IRAM geocentric apparent po sitions of the peak of 
the continuum from A ltenhoff et al. I (ll999h are given in Table |4] 
The first set of residuals shown are from Orbit I, the 2693- 
observation orbit mentionned above. Although these residuals 
are systematically negative, the important point is that those in 
declination are only slightly negative, suggesting that the peak of 
IRAM continuum emission is substantially closer to the comet's 
centre of mass than the optical astrometric observations are dur- 
ing this time. The second set of residuals are from an orbit so- 
lution (Orbit II in Table[3j that also includes these IRAM obser- 
vations. The consistency of the IRAM data suggests it may be 
appropriate to weight them more heavily in the solution, specif- 
ically, by a factor of five. At the same time, it is reasonable to 
accept in the solution only residuals in either coordinate that 
are no more than 1.5". The resulting orbital elements (Orbit 
III) are given in Table [3] These satisfy 2249 observations with 
mean residual 0.71". The residuals of the 1993 and 2005 ob- 
servations are changed from those of the previous solution by no 
more than 0. 1". The IRAM residuals from Orbit III are now very 
commendable. 

We conclude that there is no need to invoke the existence of 
nongravitational forces acting on this comet. The "original" and 
"future" barycentric values of 1/a are -i-0.003801 and -1-0.00557 1 
AU"', respectively, corresponding to periods of 4267 and 2405 
years. 



4. Modelling observations 

In order to press on with the data analysis, we have developed 
a model able to simulate both interferometric and ON-OFF ob- 
servations of a synthetic coma. In a first step, the brightness dis- 



tribution of molecular lines in the plane of the sky is computed 
using a radiative transfer model. This brightness distribution is 
used in the second step to compute synthetic ON-OFF line pro- 
files and to simulate velocity-resolved interferometric observa- 
tions. 



4.1. Radiative transfer model 

Computing the brightness distribution of a molecular rotational 
line in the plane of the sky is a classical radiative transfer prob- 
lem that can be easily solved out by numerical integration pro- 
vided we know: 1) the population distribution of the involved 
rotational levels; 2) the spatial and velocity distribution of the 
studied species; 3) the velocity dispersion, assumed here to be 
thermal at local temperature; and 4) the background emission, 
taken to be the 2.7 K cosmic background. It also requires an 
appropriate choice of number of cells and of cell size, both spa- 
tially and spectrally, in order to get sufficient accuracy. Detailed 
formulae of the cells contribution in the case of an optically thin 
coma can be found in Bockelee-Morvan et al. {in preparation). 
The code used in this paper is more general as it takes into ac- 
count opacity effects: in addition to the proper emission of each 
cell (at given Doppler velocity interval) it computes the cell at- 
tenuation of the signal emitted from background cells along the 
line of sight. The output of the programme is a spectral cube 
which provides the brightness distribution Fi(x,y) in the plane 
of the sky ix,y) for a number of velocity channels /. The com- 
puted map is 256x256 pixels wide, with pixels of 0.5" x 0.5", 
and velocity channels of 0.10 km s"'. 

The level populations in th e cells were computed accord- 
ing to the excitation model of iBiver et alj ( IT999.) . The model 
takes into account collisions with water molecules and electrons 
which control excitation in the inner coma: we used a total cross- 
section for de-excitation by neutral collisions of 1 -2x10"'^ cm^, 
and a x„e scaling parameter equal to 1 for the electron density. 
Radiati ve excitation by t he Sun radiation is taken into account 
for CS (iBiver et al.lll999h . Radiative excit ation is not con s idered 
for H2S and SO molecules. As detailed in Crovisier et al.l(ll991h 
and Biv er et a l. (2006), given their short photodissociation life- 
time, they do not undergo significant radiative excitation before 
being dissociated. This is particularly relevant in March 1997 
for comet Hale-Bopp, the collisional region of which was large 
at this time. The excitation model takes also into account radi- 
ation trapping, which main effect is to increase the size of the 
region where molecules are at thermal equilibrium. 

Because the field of view of the Plateau de Bure observa- 
tions was essentially sampling molecules with rotational level 
populations at thermal equilibrium, the temperature of the gas 
is an import ant parameter of the mode l. Both radio (single-dish) 

■ I 



(Bi ver et all 120021) and infrared (e.g.. iMagee-Sauer et al.lll999i: 
.DiSanti et al. 200 li) observations of comet Hale-Bopp in March 
1997 suggested gas temperatures ~ 120 K at radial distances of 
~ 10000 km from the nucleus. Evidence for a temperature in- 
crease with distance to nucleus is reported from long-sht spec- 
troscopy in the infrared (iMagee-Sauer et aO 19991: iDiSanti et al.l 
2001). Such an increase is predicted in coma hydrodynamic sim- 
ulations as a result of photolytic heating (Combietal. 199^. 
Therefore, we investigated two temperature laws in our model: 

- a constant temperature: T = 120 K, derived from the helio- 
centric dependenc e of T in comet Hale-Bopp published by 
IBiver et al] ( |2002|) . 

- a variable temperature: T — 20 K for r < 100 km, T — 
20 + 50 log {!^) for 100 km < r < 10000 km and T = 120 
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Table 3. J2000.0 orbital elements of C/1995 Ol (Hale-Bopp). Epoch 1997 Mar. 13.0 TT. 



Orbit 


T 


q 


Peri 


Node 


Incli. 


e 


1/a 




1997 Apr. (TT) 


(AU) 


n 


n 


n 




(AU-') 


I 


1.13588 


0.9141659 


130.58732 


282.47062 


89.43001 


0.9951316 


+0.0053255 


II 


1.13606 


0.9141660 


130.58734 


282.47061 


89.43002 


0.9951315 


+0.0053256 


III 


1.13623 


0.9141657 


130.58735 


282.47061 


89.43003 


0.9951313 


+0.0053258 




(± 0.00003)" 


(± 0.0000002)" 


(± 0.00001)" 


(± 0.00001)" 


(± 0.00001)" 


(± 0.0000002)" 


(± 0.0000002)" 



Errorbars for orbit III. 



Table 4. Astrometry of C/1995 Ol (Hale-Bopp) from continuum observations at Plateau de Bure interferometer and residuals with 
respect to orbits in Table [3] 



1997 Mar. UT 


RA (J2000.0) 


Dec (J2000.0) 


(O- 


-C) Orbit I 


(O- 


-C) Orbit II 


(O-C) Orbit III 










((5RA,(5Dec) 


(SRA,SDec) 


(SRA,6Dec) 


day 


h:min:s 


h:min:s 






(",") 




(",") 


(",") 


09 


05:00:00 


22:14:44.864 


+39:22:44:07 


(- 


1.4,-0.6) 


(- 


-0.9,-0.4) 


(-0.5,-0.1) 


09 


08:00:00 


22:15:40.688 


+39:27:56.90 


(- 


1.2-0.6) 


(- 


-0.8,-0.4) 


(-0.3,-0.1) 


09 


10:00:00 


22:16:18.031 


+39:31:24.84 


(- 


1.2,-0.6) 


(- 


-0.7,-0.4) 


(-0.2,-0.1) 


09 


12:00:00 


22:16:55.459 


+39:34:52.18 


(- 


1.4,-0.7) 


(- 


-1.0,-0.5) 


(-0.5,-0.3) 


11 


06:00:00 


22:30:26.714 


+40:45:22.83 


(- 


1.3,-0.6) 


(- 


-0.8,-0.5) 


(-0.3,-0.2) 


11 


06:00:00 


22:30:26.724 


+40:45:23.05 


(- 


1.2,-0.4) 


(- 


-0.7,-0.2) 


(-0.2, 0.0) 


13 


08:00:00 


22:47:34.109 


+42:03:01.56 


(- 


1.3,-0.5) 


(- 


-0.7,-0.4) 


(-0.2,-0.2) 


13 


15:00:00 


22:50:03.197 


+42:13:15.05 


(- 


1.3,-0.5) 


(- 


-0.7,-0.3) 


(-0.2,-0.2) 


16 


16:00:00 


23:17:10.406 


+43:48:28.38 


(- 


1.1,-0.4) 


(- 


-0.5,-0.4) 


( 0.0,-0.2) 


16 


16:00:00 


23:17:10.407 


+43:48:28.42 


(- 


1.1,-0.4) 


(- 


-0.5,-0.3) 


(+0.1,-0.2) 



" Continuum geocentric positions from lAltenhoff et 5] (Il999i) converted from apparent to J2000.0 coordinates. 



K for r > 10000 km, where r is the radial distance from the 
nucleus. 

The variable temperature law provides a temperature of 70 
K at 1000 km from the nucleus, consistent with the model of 
ICombi et al.l (11999 '). but somewhat below the rotational tem- 
peratures of 80-90 K deriv ed from infrared measurements with 
1" (1 000 km) aperture (M agee-Sauer et al.lll999t iDiSanti et al] 
l2001i) . Therefore, this temperature law allows us to investigate 
extreme temperature gradients. The population of the upper and 
lower levels of the CS, H2S and SO observed transitions are 
shown in Figs. [TOtfTTI for the two temperature laws. As the dis- 
tance from the nucleus increases, the rotational population dis- 
tributions of the molecules evolve from thermal equilibrium to 
a cold fluorescence equilibrium, because of declining collision 
rates. This explains the increase in the outer coma of the popula- 
tion of the levels involved in the observed rotational transitions. 
With this temperature law, the modelled ON-OFF observations 
(10000-25000 km field of view radius) mostly probe collision- 
ally excited molecules at T = 120 K, while the interferometric 
beam (1000 km equivalent radius) is sensitive to colder regions 
at r ~ 70 K. Colder gas in the inner coma results in larger opac- 
ity effects, as it will be shown in Sect. |5] 

The radiative transfer code can accommodate any 3-D den- 
sity and velocity distribution as inputs. For the purpose of 
Sect.|5] namely the analysis of the radial distribution of CS, H2S 
and SO, the density distribution was described by the Haser for- 
mulas for parent and/or daughter species. Since the H2S coma 
does not show up strong spatial and spectral asymmetries, the 
assumption of isotropic outgassing is satisfactory for analysing 
the H2S data. This is less true for SO and CS (Sect ID, so that 
this must be kept in mind in the discussion of the results. We can 



note that the Haser formula for daughter molecules is appropri- 
ate for describing the density distribution of SO and CS in comet 
Hale-Bopp within the IRAM instrumental field of view given the 
large size (> 10^ km, Combi et al. 1999) of the collision sphere 
in March 1997 compared to the scalelengths of their presumed 
parents (see Sect. [Sj. Numerical simulations of asymmetric co- 
mas discussed in Sect l2.3.2l were performed using a combination 
of conical structures, each of them characterized by a production 
rate by unit of solid angle. 

For the gas expansion velocity, we assumed a value about 
1.03 km s ', derived from the heliocentric dependen ce of the gas 
expan sion velocity measured in comet Hale-Bopp bv lBiver et al.l 
(l2002h . Other relevant parameters in the Haser model that will be 
constrained by the observations are the photodissociation rates 
of CS, H2S, SO and SO parent, and the CS, H2S and SO produc- 
tion rates. 



4.2. ON-OFF synthetic spectra 

Synthetic ON-OFF line profiles are readily obtained from the 
spectral cube by convolving the brightness distribution obtained 
for each velocity channel with the primary beam pattern. The 
beam pattern A(x,y) was described by a 2D-gaussian which 
width at half maximum is equal to the HPBW given in Table [T] 



4.3. Synthetic interferometric observations 

The simulation of interferometric observations is more complex, 
as it requires the Fourier transform of the brighness distribution 
taking into account the uv-coverage (see Sect. l2.3l l. 
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Table 5. Published photodissociation rates at 1 AU from the Sun. 



Table 6. Observed and modelled fluxes (in Jy km s ) and flux 
ratios. 



Molecule 



Reference 





Z.J X lU 


-4 


Crovisier et al. 1991 




i.U X iU 


-5 


JacKson et at. tyoz 




Z.U X lU 


-5 


Diver et al. zuUj 




1.0 X 10 


-4 


Snyder et al. 2001 


CS2 


1.0 X 10 


-3 


Feldman et al. 1999 




1.7 X 10 


-3 


Jackson et al. 1986 


SO 


1.5 X 10 


-4 


Kim& A'Hearn 1991 




4.9 X 10 


-4 


Summers & Strobel 1996 




6.2 X 10 


-4 


Huebner et al. 1992 


S02 


2.1 X 10- 


-4 


Kim & A'Hearn 1991 




2.1 X 10- 


-4 


Huebner et al. 1992 




2.3 X 10 


-4 


Kumar 1982 




2.9 X 10 


-4 * 


Summers & Strobel 1996 



Observations" 



Models* 



* The solar spectrum model used pertains to solar maximum. 



For each channel flie visibilities are defined by (see e.g., 
[Thompson et~an.ll 9911) : 



%(o-) 



- r 



A(s)Fi(s) exp(-^^o- ■ s)dn, 

4n <^ 



(1) 



where cr is the baseline vector for two antennas, with coordinates 
(m,v) in the Mv-plane. s is a vector in the sky plane which coor- 
dinates are (x,y) in radian units. A is the power pattern of the 
antennas, and d£l is an element of solid angle on the sky. Unit 
of the brightness distribution F,(x, y) provided by the radiative 
transfer model is [W m"^ sr"']. 'V,(cr) is here in units of [W m"^ 
Hz '] or janskys. Equation[T]can be approximated to : 

00 U —00 



v6v 



A(x,y)Fi{x,y) 



X exp( —{ux + vy)) dxdy. 



(2) 



The calculation of the visibilities is carried out using the 
GILDAS software collection provided by IRAM. The first step 
consists in converting the brightness spectral map Fjix, y) times 
the primary beam gaussian pattern A{x, y) into GILDAS format. 
Then a GILDAS task computes the Fourier transform accord- 
ing to Eq. 12] In order to simulate data comparable to observa- 
tions, it is necessary to take into account the uv-plane sampled 
by the observations. In practice, the GILDAS task reads the uv- 
coverage in the observed data file and computes the 'Viiu, v) for 
each (m,v) point. This step provides a wv-table that can be handled 
exactly like that recorded during the observations. Modelled and 
observed maps obtained by inverse Fourier Transform can be 
compared. 

It is also possible to compare the visibilities in the Fourier 
space. Because of symmetry properties of the Fourier Transform, 
a radially symmetric brightness distribution results in visibili- 
ties 'Vi{u, v) that only depend upon the wv-radius cr - Vm^ + v^. 
Observed visibilities that have been radially averaged in the uv- 
plane can be directly compared to modelled 'V,(cr). In such a 
way, the mean (angularly averaged) radial brightness distribu- 
tion of observed species can be investigated. For this study, visi- 
bilities integrated over velocity 2; 'Vi{a-)6v [Jy km s"'] (with 6v 
equal to channel width) will be compared (Sect.|5]). Note that the 
line flux measured in ON-OFF mode corresponds to the visibil- 
ity at cr = 0. 







T = 


120 K 


T variable 






Thin 


Thick 


Thin 


Thick 


H2S 220- 


-2,1 










FsD 


18.9 ±2.8 


19.14 


17.20 


19.40 


16.60 


Pint 


3.51 ±0.55 


4.89 


3.96 


5.15 


3.60 


R 


5.38 ± 1.16 


i.yi 




3.76 


4.61 


CS J(2- 


■1) 










FsD 


15.2 ±2.3 


15.1 


14.3 


14.6 


13.7 


Pint 


1.47 ±0.24 


1.07 


1.02 


1.50 


1.29 


R 


10.3 ±2.3 


14.1 


14.0 


9.7 


10.6 


CS J(5^) 










FsD 


213 ±32 


231 


213 


227 


201 


Pint 


5.4 ± 1.1 


11.4 


8.22 


14.7 


7.30 


R 


39± 10 


20.3 


25.9 


15.4 


27.5 



Line integrated intensity in ON-OFF mode (Fsd) and at the cen- 
tre of the interferometric map (F|nt) and their ratio R = Fsd /Pint- 
Unlike in Table[T] error bars take into account a 15% uncertainty in 
the calibration. 

Model calculations of Fsd and Fi„t, assuming a constant or vari- 
able temperature in the coma (see Sect. Results labelled as 
thick (respectively thin) indicate whether opacity effects were in- 
cluded or not in the calculations. Calculations were performed with 
/3h,s = 2.5 X lO-"* .r',/3cs, = 1-7 x IQ-^ s-\ /3cs = 2 x IQ-^ s"'. 
The assumed production rates are (2h,s = 1.35 x 10^' s"', Qcs = 
2.0 X 10-** S-' (respectively 2h,s = L2"x 10-' s"', 2cs = 1.8 x 10^* 



s ) for T = 120 K (respectively, variable T) law. 



5. Analysis of the radial distribution of CS, H2S and 
SO 

Two different methods are used to analyse the radial distribu- 
tion of CS, H2S and SO. Let us define F^d as the line area of 
the ON-OFF spectra and Fi„t as the line area at maximum in- 
tensity in the interferometric maps, which observed values are 
given in Table [1] Because of the factor of ~ 10 difference in 
angular resolution between ON-OFF and interferometric obser- 
vations, the flux ratio R = /^so/^int is strongly sensitive to the 
spatial distribution. We can thus compare observed and modelled 
flux ratios to constrain model parameters. The more extended the 
brightness distribution is, the larger the flux ratio R. In the sec- 
ond method, we compare observed and modelled visibilities as 
a function of Mv-radius, as explained previously, including the 
ON-OFF data point. 

Input free parameters in the model are photodissociation and 
production rates. The production rate Q is adjusted iteratively 
by running the model several times in order to match at best the 
flux ratio and the observed intensities Fsd and Fint- Indeed, the 
intensities (particularly Fi„t) are aff'ected by small but significant 
opacity effects, so that the flux ratio depends on Q. 

5.1. Hydrogen sulfide H2S 

The photodiss ociation rate of H2S has been determined by 
ICrovisier et al.l d 199 lb using laboratory measurements (Table |5]l. 
The resulting H2S scalelength for comet Hale-Bopp in March 
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1997 is ~ 3900 km, i.e. intermediate between the angular reso- 
lution of ON-OFF and interferometric observations. Therefore, 
the flux ratio R is expected to reflect, in large part, H2S photoly- 
sis. 

The observed flux ratio R can be reproduced within 1-cr with 
the model for a photodissociation rate of H2S in the range 1 .5 - 
3x lO^^s^'.This range includes the value from Crovisi er et al.l 
( Il991h . for which computed flux densities and flux ratios R are 
given in Table |6] for the two temperature laws, together with 
measured values. Calculations in the assumption of an optically 
thin line are also given. The good agreement between model and 
observations (within l-cr) in the optically thick case shows that 
the radial extension of H2S in the coma is well explained by 
direct release of H2S from the nucleus. Opacity effects signifi- 
cantly affect the line flux in the interferometric beam: by 23% 
with the constant temperature law to 40% for T variable. The 
retrieved H2S production rates are 1.35 and 1.2 x 10^^ s"', for 
the constant and variable T laws, respectively. 

Figure [12] shows the comparison between modelled (with 
opacity effects included) and measured visibilities for the two 
temperature laws. Both models reproduce correctly the radial 
evolution of the visibilities. At large Mv-radii, the fit is better 
with the model assuming a variable temperature. Fluctuations 
apart the model fit trace spatial asymmetries in the coma. 

The ON-OFF H2S spectrum given by the model is super- 
imposed on the observed spectrum in Fig |3] The agreement is 
satisfactory, not considering the small spectral asymmetry in the 
observed spectrum which is due to non-isotropic H2S density 
and velocity distribution as discussed in Sect. |2] 

5.2. Carbon sulfide CS 

Ultraviolet observations of CS performed with the International 
Ultraviolet Explorer (lUE) have shown that CS2 is the likely par- 
ent of CS dJackson et al.lll982lll986l) . More recently, a compar- 
ison between unidentified lines in comet 122P/de Vico spectra 
and jet cooled laboratory CS2 spe ctra suggests that CS 2 is indeed 
present in cometary atmospheres ^Jackson et al.'2004). The CS2 
lifetime at 1 AU is estimated betw een 590 s (Jackson et al. 1986) 
and 1000 s jFeldman et al.ll 19991) from observations and labora- 
tory measurements (i.e., photodissociation rate yScs, between 1 
and 1.7 x 10"^ s"'). 

The ph otodissociation rate of CS is debated. From lUE ob- 
servations, IJacksonjeOiP (fT982) deduced yScs = 1 x 10"^ s 
From CS /(2-1) interferometric obser vations of comet H ale- 
Bopp performed with the BIMA array, ISnvder et al] (1200 ll) in- 
ferred yScs = 1 X 10 s i.e. a value ten times higher On 
the other hand, single-dish millimetric observations of comets 
at small heliocen tric distances suggested /3cs ~ 2 x 10"^ s"' 
dBiver et al ] l200l . 

The field of view radius of our CS ON-OFF observations is 
10000 to 25000 km, depending on the line (TablelTJ. Therefore, 
the calculated flux ratio Fsd and Fim should be strongly depen- 
dent upon the assumed CS scalelength, when exploring values 
between 10000 and 100000 km. On the other hand, the interfer- 
ometric beams (750 to 2000 km radius, for the J(5-4-) and 7(2-1 ) 
lines, respectively) are comparable to the CS2 scalelength, which 
makes the determination of the CS2 lifetime more difficult, given 
large uncertainties in CS lifetime. 

Figure[T3]shows the dependence of the flux ratio of the 7(2- 
1 ) (R21 ) and J(5-4-) {R54) lines with the photodissociation rate of 
CS, assuming a constant or a variable temperature in the coma. 
The CS2 photodissociation rate was taken equal to y6cs2 = 1-7 
X lO"-' s"' (value at r/, - 1 AU). The observed ratios can be 



reproduced within 2-cr, with /3cs in the range 2-7 x 10"^ s"' 
(case of constant T) or with /3cs between 0.7 and 6 x 10"^ s"' 
(case of variable T) . The same co mputations were made with 
/3cs. = 1 X 10"^ s"' (iFeldman et al.|[l999h in order to investigate 
the sensitivity of the results to this parameter. A satisfactory fit 
within 1-0" is obtained for both lines {R21 - 12 and/?54 = 29) with 
jScs in the range 2-5 x 10"^ s"' and a variable T. In summary, 
our data are consistent with a CS photodissociation rate of 1- 
5 X 10"'' s"^ On the basis of the 7(2-1) data alone, the value 
yScs - 2-3 X 10"^ s"' is preferred. Models assuming a variable 
temperature in the coma provide the best results. 

Table |6] shows that opacity effects strongly affect the inten- 
sity of the 7(5-4) line in the interferometric beam, justifying our 
radiative transfer model. They are less significant for 7(2-1), and 
for r = 120 K than for variable T. These trends can be easily ex- 
plained: 1) the 7(5^) is intrinsically stronger and observed with 
a smaller field of view; 2) with the adopted variable tempera- 
ture law, the lower rotational levels of the 7(2-1) and 7(5^) 
transitions (namely 7=1 and 7 = 4) are more significantly pop- 
ulated in the inner coma (see Fig. [TOl i. hence more significant 
self-absorption effects. The retrieved CS production rate is ~ 2 
X 10-** S-' withyScs = 1-2 X 10"^ s-\ 

In Fig. [14] we show the radial evolution of the visibilities 
for different model parameters, assuming a variable temperature 
in the coma. Models with /5cs - 2 x 10"^ s"' and Pcs^ equal 
to 1 X 10""^ s"' and 1.7 x lO""^ s"' reproduce correctly the ob- 
served visibilities, as expected from the study of the fluxes ratio. 
Visibilities computed with /3cs = 1 x 10"'* s"' do not match the 
data. 

Th e high value Pes = 1 x 10 ''s ' inferred by Snyd er et al.l 
(1200 ll) from comet Hale-Bopp observations of CS 7(2-1) with 
the BIMA arr ay is ruled ou t by ou r analysis of Plateau de Bure 
observations. ISnyder et al.l (1200 lb used a different approach. 
They analysed CS radial profiles deduced from the interferomet- 
ric maps, but did not consider the autocorrelation measurements 
that they acquired nor the radial evolution of the visibilities in 
the Mv-plane. However, for yScs = 1 x 10~*s~', their autocor- 
relation data yield systematically higher (by a factor 2.5-3) CS 
production rates than the cross-correlation data. We computed 
that, if yScs is decreased to ~ 3 x 10"^ s"', the production rates 
derived by the two sets of data can be reconciled. 

5.3. Sulfur monoxide SO 

As pointed out in Sect. 12.3.21 there is clear evidence that the 
SO radial brightness distribution is much more extended than 
that of H2S. The flux ratio for the SO line is = 38 + 10, 
taking into account calibration uncertainties, to be compared 
to R = 5.4 ± 1.2 for the H2S fine (Table |6ll. SO is a radical 
and a photodi ssociation product of SO2, a lso observed in comet 
Hale-Bopp (B ockelee-Morvan et al.l l2000^. Published SO2 pho- 
todissociation rates are in general consistent with each other: 
^so, = 2.1 - 2.3 X lO""* s-' for quiet sun and 2.9 x 10"^ ^-i foj. 
act ive sun ( see Ta ble[5]l. A significantly lower value was derived 
bv lKumad ( |1982|) . but only one channel of the SO2 photodisso- 
ciation process was considered. We will not consider this value 
in this work. 

In contrast, published values of Pso differ by up to a factor of 
4 (Table|5]l. Because SO and SO2 photodissociative scalelengths 
are intermediate between the primary and interferometric fields 
of view, they can be somewhat constrained by the observations. 
Model calculations were performed for many couples (ySso, ySso, ) 
in the range of published values. Isocontours of the flux ratio R 
as a function of Pso and Pso^ are presented in Fig. [15] 
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In the assumption that SO is produced solely by the 
photodissociation of SO2, the observations are consistent 
(within 1 -1.5-cr) with the pub lished value of ySso - 1.5 x 
10 "^s ' (iKim & A' Heard 1 199 lb . but ar e not compatib l e with 
values as high as 5-6 x lO^'^s"' jHuebner et all 119921; 
[Summers & Strobell Il996l) . The same conclusion is obtained 
from the study of the visibility curve (Fig. [T6] l. Direct release 
of SO from the nucleus is fully excluded (Fig.fTSIl. 

Published ySso rates are all based on the same laboratory data 
(Phillips, 1981). Discrepancies may be ascribed to the complex- 
ity of the laboratory SO absorption spectrum (which shows a 
series of strong narrow lines over a broad continuum) and the 
fact that the spectrum is only available in a graphic version. 
From our own analysis of the spectrum, we put a firm lower 
limit j6so > 1.5 x 10"^ s"' using the lower envelope to the SO 
spectrum. Considering in addition the narrow lines, we estimate 
jSso = 3.2 X 10""* s"', but this estimation is still uncertain due to 
the lack of a digital version of the spectrum. 

How this result could be reconciled with our measurements? 
First, because SO emission is strongly anisotropic, can we ex- 
clude a flaw in our analysis due to simple approximations made 
on the coma structure? Calculations performed for an anisotropic 
coma which reproduces the main characteristics of the SO map 
(Sect l2.3.2b show that our analysis should not be significantly 
affected by the presence of jets in the coma. In addition, a sepa- 
rate reduction of the blue and red velocity parts of the SO spec- 
trum shows that both parts display the same intensity flux ratio R 
(equal to that measured for the whole emission), indicating that 
jet structures do not affect the overall evolution of the visibilities 
with baseline length. 

A /3so value higher than that which fits the observations im- 
plies that some process makes the SO brightness radial distribu- 
tion more extended than expected from SO2 photodissociation. 
We examined several possibilities. Among those that can be ex- 
cluded are: 1) a ^6502 value lower than published values; shield- 
ing of SO2 photodissociation by water is insufficient (the inner 
coma where UV shielding by water is significant may be eval- 
uated to r < 1000 km); 2) a kinetic temperature that is much 
higher ( ~ 200 K) in the inner coma observed by the inter- 
ferometric beam than in the outer coma observed by the pri- 
mary beam, that would be in contradiction with measurements 
dMapee-Sauer et al.L 119991; iDiSanti et all I2OOII) and hydrody- 
namical models dCombi et al .L I1999I) . This leaves the possibility 
that SO2 is not the only source of SO, or that SO2 is itself pro- 
duced by some extended source in the coma such as grains or 
complex molecules. 

The SO production rate deduced from the observations 
is 2.4-2.7 X lO^** s-i for ^650 = 1-5 x lO-'^s"' and 
Pso. = 2.1 X 10~'*s~ '. Base d on the same ON-OFF data, 
iBockelee-Morvan et al. I (120001) inferred a higher value for 
the SO production rate due to different assumptions on the 
beam efficiency and kinetic temperature. SO2 was observed 
in comet Hale-Bopp on several dates in March and April 
1997 using the IRAM Plateau de Bure and 30-m antennas 
(iBockele e-Morvan et al."2000). The measurements yield a pro- 
duction ratio 2so/2sOt ~ 2 in March 1997. The discrepancy in- 
creases when a higher fis p value is adopted (e.g., Qso/Q so, ~ 6 
for/3so = 6.2x 10'^ s-^). lBockelee-Morvan et alJdloOOh suggest 
that Qso^ is possibly underestimated in their calculations as they 
assumed LTE for the SO2 rotational population distribution, and 
neglect electronic and vibrational excitation. Another possible 
explanation is the presence of another source of SO in the coma. 
As discussed above, this would allow a more satisfactory inter- 



pretation of the observed SO radial distribution, given available 
estimations of /3so. 

The synthetic SO ON-OFF spectrum computed with our 
isotropic model is plotted in Fig. [3] Though it does not match 
in details the observed asymmetric SO spectrum, the line width 
and wings are correctly reproduced. This indicates that, as as- 
sumed in the model, the SO radicals which are within the field 
of view, though daughter products, do indeed partake in the ki- 
netics of the gas flow. 

6. Summary 

Observations carried out in cross- and autocorrelation modes 
with the IRAM Plateau de Bure interferometer were analysed to 
study the spatial distribution of sulfur-bearing molecules in the 
coma of comet Hale-Bopp using a radiative transfer code. The 
radial distribution of H2S is in agreement with a direct release 
from the nucleus and is well reproduced by our model using the 
H2S photodissociation rate published in the literature, thereby 
validating our approach. 

On this basis, we studied the radial distribution of CS emis- 
sion. This radical is believed to be produced by the photolysis of 
the short-lived CS2 molecule. The angular resolution of the ob- 
servations prevented us to obtain observational constraints on the 
CS2 photodissociation rate. The photodissociation rate of CS is 
not precisely determined because no satisfactory measurements 
can be done in the laboratory on this radical. Our data are consis- 
tent with j6cs in the range 1 - 5 x 10"^ s"', which is in agreement 
with the value of 2 x 10"^ s"' proposed by Biver et al. (2003) 
but rules out yScs = 1 x lO^'^s"' derived bv lSnvder et al.1 020011) 
from interferometric observations with the BIMA array. With 
/3cs ~ 3 X 10"^ s the CS production rates derived from auto- 
correlation and cross-correlation BIMA data are reconciled. 

SO was first detected in co met H ale-Bopp, as we ll as its 
presumed parent product SO2 (iBockelee-Morvan et all l2000l) . 
Values of /3s o available in the litera ture differ by a factor of ~4 
between the lowest (1.5 x 10'"* s~', iKim & AH earn 1991) and 
the highest value (6.2 x 10"* s '.'Hue bner et al.fl992) . From our 
own analysis based on the same laboratory absorption spectrum, 
we estimate /3s o ~ 3.2 x 10""* s"'. In contrast, assuming that 
SO2 is the parent of SO, the SO interferometric data suggest /^so 
= 1.5 X 10~^ s"'. The discrepancy between measured and com- 
puted values of the SO photodissociation rate may indicate that 
SO2 is not the sole parent of SO, or that SO2 is itself created by 
some extended source in the coma. Note that the sulfur chem- 
istry in comets is far from being fully u nderstood: the origins of 
S2 a nd NS radicals are still unknown (iBockelee-Morvan et all 
'2004). In the future, interferometric mapping of SO and SO2 
with the Atacama Large Millimeter and submillimeter Array 
(ALMA) will provide relevant informati ons about the SO cre- 
ation process in cometary atmospheres ( Biverll2005h . Further 
laboratory experiments on SO absorption in the UV are also 
strongly encouraged. 

The data show clear evidence for anisotropic spatial distribu- 
tions in the coma: (1) ON-OFF spectra are asymmetric; (2) the 
visibilities as a function of Mv-radius do not follow the expected 
curves in the case of a 1-D outflow ; (3) maps of CS and SO 
line brightness distributions do not peak at the position of the 
continuum maps obtained simultaneously which should be close 
to the nucleus position; (4) maps show strong and specific spec- 
tral variations in the coma. On the basis of model calculations, 
we show that SO and CS present a wide jet-like structure that 
may correspond to the gaseous counterpart of the strong high 
latitude dust jet detected at optical wavelengths. This structure is 
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not observed for H2S. All three molecules show day/night spa- 
tial asymmetries. Spectral temporal variations related to nucleus 
rotation are only observed for CS. The strong CO rotating jet 
observed at Plateau de Bure interferometer has no H2S, CS and 
SO counterpart. These differences in the outgassing patterns of 
different volatiles in the Hale-Bopp coma may suggest that the 
nucleus composition is inhomogeneous (though all molecules 
were not observed on the same day and we cannot exclude time 
variations of the outgassing). This analysis is consistent with 
the results of some previous works: (1) the offsets observed be- 
tween molecular (HNC, DCN, HDO) and continuum emission 
on interferometric maps made with the Owens Valley Radio 
Observatory (OVR O) were interprete d as due to jets enriched 
in these molecules (iBlake et al.Lll999b ; (2) Monte Carlo simula- 
tions suggest that the different molecular jets observed at optical 
wavelengths (OH, CN, C2) origi nate in several activ e areas of 
different chemical compositions dLederer et al.L 120021) . An ex- 
planation would be that the huge nucleus of comet Hale-Bopp is 
actually made of several pieces that have not undergone the same 
formation and/or evolution processes. This would be in contrast 
with comet 73P/Schwassmann-Wachmann 3 whose fragments 
were found simila r in chemical properties, suggesting an homo- 
geneous nucleus (iDello Russo et al.ll2007l) . 

The results obtained in this paper demonstrate the benefit 
and uniqueness of millimeter interferometry to investigate nu- 
cleus outgassing properties and physical processes acting in the 
coma. Last and not least, by solving the contentious issue regard- 
ing the role of nongravitational forces in comet Hale-Bopp orbit, 
we have shown that this technique is also powerful for astrome- 
try purposes. ALMA, with its unprecedented sensitivity, angular 
resolution, and instantaneous Mv-coverage, will be a very power- 
ful tool in the future for cometary science. 
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Fig, 1, a-b: ON-OFF spectra of CS 7(5^) (244.9 GHz) (a) and 
CS 7(2-1) (98.0 GHz) (b) lines recorded on March 12, 1997. 
The spectra, in brightness temperature scale, are shifted verti- 
cally according to observation UT time (given in the top right 
corners). Horizontal lines correspond to zero intensity. The in- 
tegration time is 2 min (CS 7(5-4) at 6.40 h and 14.04 h and 
CS 7(2-1) at 8.69 h and 14.04 h) or 4 min. The velocity scale is 
with respect to the comet rest velocity, c: Line area as a fonction 
of UT time. Black and grey symbols are for 7(5^) and 7(2- 
1) lines, respectively. The black dashed (respectively grey dot- 
ted) line shows mean value for 7(5^) (respectively 7(2-1)) line. 
d: Line velocity shift as a fonction of UT time. Black and grey 
symbols are for 7(5-4) and 7(2-1) lines, respectively. The black 
dashed (respectively grey dotted) line shows mean Av value for 
7(5-4-) (respectively 7(2-1)) Une. 
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Fig. 2. a-b: ON-OFF spectra of H2S 220 - 2n (216.7 GHz) (a) 
and SO NjiSe-^s) (219.9 GHz) (b) hues recorded on March 13, 
1997. The spectra, in brightness temperature scale, are shifted 
vertically according to observation UT time (given in the top 
right corners). Horizontal lines correspond to zero intensity. The 
integration time is 4 min for all spectra, but those at 4.32 h and 
14.62 h UT for which it is 2 min. The velocity scale is with re- 
spect to the comet rest velocity, c: Line area as a fonction of 
UT time. Black and grey symbols are for H2S and SO, respec- 
tively. The black dashed (respectively grey dotted) line shows 
mean value for H2S (respectively SO) line, d: Line velocity shift 
as a fonction of UT time. Black and grey symbols are for H2S 
and SO, respectively. The black dashed (respectively grey dot- 
ted) Une shows mean Av value for H2S (respectively SO) Une. 
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Fig. 8. Same figure as Fig.|6]for CS 7(2-1) (top) and 7(5-4) (bot- 
tom). The integration time is 4.7 and 5.9 h, respectively. At pos- 
itive Dec offsets (and especially at positive RA), the blueshifted 
side of the lines is more intense; spectra are globally symmet- 
ric at negative Dec offsets. In addition, at zero and negative Dec 
offsets, 7(5-4) spectra present a narrow component peaking at 
V ~ 0. Note that the lack of flux in the Dec=-2" row of CS 7(5- 
4) spectral map may be due to the lack of short spacings in the 
Mv-plane. 
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Fig. 9. Synthetic spectral maps (left panel), ON-OFF spectra 
(right upper panel) and positions of the brightness maximum 
with respect to the nucleus (right lower panel) for CS /(5^) 
(top) and 7(2-1) (bottom) lines as obtained with a model includ- 
ing a jet like structure southward and a day/night asymmetry in 
the coma. The jet is 50° wide, and contains about 25 % the total 
amount of CS. The CS production rate is 5 times stronger in the 
day side of the nucleus. The arrow indicates the projection of the 
Sun direction in the plane of the sky. The phase angle is 45°. 
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Fig. 10. Population of the rotational levels of CS related to ob- 
served transitions: J = 1, 2, 4 and 5 as computed with the excita- 
tion model. The upper panel corresponds to a constant tempera- 
ture, the lower to a variable temperature. 
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Fig. 11. Population of the H2S and SO rotational levels related to 
observed transitions as computed with the excitation model. The 
upper panel corresponds to a constant temperature, the lower to 
a variable temperature. 
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Fig. 12. Amplitude of the visibility of H2S 217 GHz line ver- 
sus Mv-radius. Circles with error bars: observations. Dotted Unes: 
model assuming constant T. Thick lines: model assuming vari- 
able T. In both cases, calculations were performed with ;8h2S = 
2.5 X 10-4 s-i. 
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Fig. 14. Amplitude of the visibility as a function of wv-radius 
for CS /(2-1) (upper panel) and CS /(5-4) (lower panel). 
Measurements are shown by dots with error bars. Displayed 
modelled visibilities are for a variable temperature law and pa- 
rameter sets: (/3cs = 2 x 10"^ s ',j8cs2 = 1-7 x 10"^ s"') in solid 
line, (jScs = 2 x 10"^ s"',j8cs2 = 1 X 10"^ s"') in dashed line and 
(jScs = 1 X 10-'' s-i,jecs2 = 1-7 X 10-3 s-i) in dotted line. 




Fig. 13. Modelled flux ratios /? as a function of CS photodissoci- 

ation rate for CS /(2-1) (/?2i, top) and /(5-4) {R^^, bottom) hues 
foryScs, = 1.7 X IQ-^ s"' (value at r/, = 1 AU). The dotted (re- 
spectively dashed) lines show the 1-cr (resp. 2-cr) range around 
the observed value (continuous Une). 
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Fig. 15. The flux ratio R - ^ as a function of f3so and /3s 02- 
Upper panel: for constant temperature T = 120 K; lower panel: 
for the variable temperature law. Contour values are given in 
the bottom of the figures. Circles correspond to (J3so, PsOt) 
ues published in the literature (Table |5]). The region where the 
computed R is in agreement within 1-cr with the observations 
(/? = 38 ± 10) is filled in grey. 
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Fig. 16. Amplitude of visibility of SO 220 GHz line versus uv- 
radius. Observations: circles with error bars. Continuous, dot- 
ted and dashed lines correspond to model calculations assum- 
ing that SO is created by the photodissociation of SO2 with 
Pso of 1.5, 4.9 and 6.2 x lO^'^s"', respectively. The dotted- 
dashed line tests the case of a pure nuclear origin for SO, with 
Pso = 1.5 X 10""* s"' . Calculations correspond to variable T law. 
Similar results are obtained with a constant T - 120 K. 



